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Abstract 

We investigate the effect of K~ condensed matter on bulk viscosity and r-mode instability in 
neutron stars. The bulk viscosity coefficient due to the non-leptonic process n p + K~ is studied 
here. In this connection, equations of state are constructed within the framework of relativistic 
field theoretical models where nucleon-nucleon and kaon-nucleon interactions are mediated by the 
exchange of scalar and vector mesons. We find that the bulk viscosity coefficient due to the non- 
leptonic weak process in the condensate is suppressed by several orders of magnitude. Consequently, 
kaon bulk viscosity may not damp the r-mode instability in neutron stars. 

PACS numbers: 97.60.Jd, 26.60,+c, 04.40.Dg 
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I. INTRODUCTION 



Astrophysical compact objects such as neutron stars and black holes, are candidates of 
gravitational waves. Neutron stars pulsate in different modes due to its fluid perturbation. 
The Coriolis restored inertial r-modes of rapidly rotating neutron stars may be sources of 
detectable gravitational waves. Gravitational radiation drives the r-modes unstable due to 
Chandrasekhar-Friedman-Schutz mechanism l[ [2], [3, 0, 0, (], Q, [9] . On the other hand, this 
instability may be suppressed in different ways. One plausible explanation for the damping 



of the instability is the large bul 



hyperons 10, 
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w. 
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if hyperons are superfluid 



i viscosity coefficient due to non-leptonic processes involving 



laliz 



6 [.However, it may not be an efficient damping mechanism 
, |l8( . Another possibility of damping the r-mode instab ility 
is the mutual friction between inter-penetrating neutron and proton superfluids 18|, [19]. 
This shows that superfluidity of neutron star matter plays a significant role on the damping 
of the r-modes in neutron stars. 

Besides hyperons, other novel phases with large strangeness fraction such as quark matter 
and Bose-Einstein condensed matter of K~ mesons, may appear in neutron star interior. 



Already the bulk viscosity of unpaired guar 
ity were investigated by various groups 



i matter and its influence on the r-mode instabil- 
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231 ] . Recently, bulk viscosity coefficients 



due to weak processes in quark-hadron mixed phase 
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quark matter 
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and color superconducting 



have been studied by several groups. However, there are no 
calculations of bulk viscosity in Bose-Einstein condensed matter so far. 

Kaplan and Nelson 29] first demonstrated that Bose-Einstein condensation of i^~mesons 
in dense hadronic matter was a possibility. This calculation was performed using a chiral 
SU (3) l x SU (3) r Lagrangian. It was noted that the strongly attractive i^~-nucleon in- 
teraction lowered the effective mass and in-medium energy of K~ mesons in dense mat- 
ter. The s-wave condensation sets in when the in-medium energy of negatively charged 
kaons equals to its chemical potential. Kaon condensation in neutron star interior was 



inve stig ated extensi vely in the chiral 30, 
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and the meson exchange models 
421 ] . The onset of K~ condensation in neutron star mat- 



ter depends on the nuclear equation of state as well as on the depth of antikaon optical 
potential. Earlier calculations showed that kaon condensation might set in at around two 
times normal nuclear matter density 42|. However, the early appearance of hyperons may 
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delay K~ condensation to higher density or vice versa. 

With the onset of K~ condensation, the neutron star matter is in chemical equilibrium 
through the weak process n p + K~ . The density perturbation associated with the r- 
modes drives the system out of chemical equilibrium. The weak non-leptonic processes may 
reestablish the equilibrium. Therefore, it motivates us to investigate the bulk viscosity due 
to this non-leptonic process in kaon condensed matter and its influence on neutron star 
r-modes. The paper is organised in the following way. In section 2, we describe the field 
theoretical models of strong interactions, different phases of dense matter, the bulk viscosity 
coefficient and the corresponding time scale. Results of our calculation are explained in 
section 3. Section 4 gives summary and conclusions. 



II. FORMALISM 



Here we consider a first order phase transition from hadronic to K~ condensed matter 
in neutron stars. Both the pure hadronic and K~ condensed matter are described within 
the framework of relativistic field theoretical models. The constituents of matter in both 
phases are neutrons (n), protons (p), electrons and muons. Neutrons and protons embedded 
in the K~ condensate behave dynamically differently than those of the hadronic phase, 
t was attributed to different mean fields which nucleons experienced in the pure phases 
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35 



421 ]. Local charge neutrality and beta-equilibrium conditions are satisfied in both 
phases. The baryon-baryon interaction is mediated by the exchange of scalar and vector 
mesons. Therefore the Lagrangian density for the pure hadronic phase is given by 

£-B = 4>B (nV*^ ~ m B + 9<tB(? ~ QuiBl^ ~ QpBl^B ' f^) IpB 

B=n,p 

+ l - (d^a - mV) - U{a) 

-\^ v + \ml^ - \ Pilv ■ pT + \m\ Pli ■ p» . (1) 
Here ipB denotes the Dirac bispinor for baryons B with vacuum mass m# and the isospin 



operator is t#. The scalar self-interaction term 



43| is 



U{a) = l -g 2 a* + \g,a A . (2) 
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We perform this calculation in the mean field approximation 44j. The mean meson fields 
a, uo and po3 in the hadronic phase are given by 

2 OU ^ — > ^ Q . . 

m a° = —fo + }^ 9 ° Bn B > ( 3 ) 
B=n,p 

m ^o = 9u>Bn h B , (4) 



m p 



B=n,p 

2 



P03 = E 9 P BhBn h B . (5) 

B=n,p 



Here n B and are scalar and baryon number densities in the hadronic phase. 
The total charge density of the hadronic phase is 

Q h = 1 Bn B ~ n - ~ n » = ' ( 6 ) 

B=n,p 

where q B is the charge of baryons B in the pure hadronic phase and n e and are charge 
densities of electrons and muons respectively. The /5-equilibrium in this phase leads to 
A*n = ^p + A*e- The chemical potential for baryons B is given by 

fJ'B = { k F B + m *B ) + 9u>BV0 + hB9 P BP03 , (?) 

where effective baryon mass is m* B = m B — g aB a and isospin projection for baryons B is 
hB- 

The energy density and pressure in this phase are given by 

h ^ 2 2 , ^ 3|1 4 I 1 2 2 I 1 2 2 

£ = 2 m ° a + 3^ 2CT + 4# 3CT + 2 mujUJ ° + 2 mpp03 

B=n,p 27T J ° l=e-,p- n J ° 



3^ 2 ^ J {k 2 +m f y/ 2 + 3 i= ^ ^J (fc2 +m 2)l/2 



(9) 



The baryon-baryon interaction in the pure kaon condensed phase has the same form as 
given by Eq. (1). In this case, we adopt a relativistic field theoretical approach for the 



description of (anti)kaon-baryon interaction 3J, |35|, |42|. Besides the exchange of a, uo and 



p mesons, ( ant i ) kaon- ( ant i) kaon interaction is also mediated by two strange mesons, scalar 



meson /o(975) (denoted hereafter as a*) and the vector meson 0(1020). As nucleons do not 
couple with strange mesons, g a *N = g<f>N = 0. The Lagrangian density for (anti)kaons in the 
minimal coupling scheme is, 

C K = D^KD^K - m%KK , (10) 

where the covariant derivative is = + ig^K^n + W<pK<f>n + W p k^k ' P p an d the effective 
mass of (anti)kaons is m* K = uik — g a K<y — g^R^*- The mean meson fields in the condensed 
phase are denoted by a, a*, u , <p and p 03 . The in-medium energy of K~ mesons for s-wave 
{k = 0) condensation is given by 

uk- = Pk- =m K - 9u>kU - g<p K <f>o + hK-9 P Kpm , (11) 

where is the chemical potential of K~ mesons and the isospin projection is I^k- — ~ 1/2. 
The mean fields in the K~ condensed phase are 



?~+ $Z 9aBnz ,S + g aK n K - , (12) 



da 

B=n,p 



2 * 

m„*a 



9a*Bnz' S + g a * K n K - , (13) 

B=n,p 

mlu = ^ 9u;BnB - g^Kn K - , (14) 

B=n,p 

™l<Po = ^2 9<t> Bn B ~ 9<$>Kn K - , (15) 

B=n,p 

m lP03 = ^2 9f>BhBTiB + 9 P KhK-n K - . (16) 
B=n,p 

The scalar and baryon number densities in the condensed phase are given by 



2vr 2 J { k 2 + m *K 2 y/2 
nf = (2J B + 1)|§, (18) 

where Fermi momentum is kp B and spin is Jb- The effective baryon mass {m B K ) in the 
condensate has the same form as in the hadronic phase. In the absence of source terms for 
K~ mesons in Eq. (12)-(16), we retain the mean fields in the hadronic phase. The scalar 



density of K mesons in the condensate is given by [35] 

n K - = 2 \ uj k - + guKOJo + gj>K<f>o + ]^9 P KPo^j KK = 2m* K KK . (19) 



The total charge density in the condensed phase is 

Q K = ^ Bn B -n- K -n e -n^ = $ (20) 

B=n,p 

With the onset of K~ condensation, we have 

n ^p + K~ . (21) 
The requirement of chemical equilibrium yields 

fj£-tf = HK-=He • (22) 

Here the baryon chemical potential has the same form as given by Eq. (7), but the mean 
fields used here are obtained in the presence of the condensate. 

The total energy density in the kaon condensed phase is comprised of baryons, leptons 
and antikaons 

K 2 2 i 3 i 4 i 2 ^2 i 2 2 i ^~ 2 j,2 i 2 2 

£ = 2 m - a +3^ 20r +2 r? V or + 2 mujUJ ° 2 m ^° 2 mpPm 

+m* K n K - . (23) 

The last term denotes the contribution of the K~ condensate. The pressure in this phase is 
1111 111 

t->K L 2 2 L 3 1 4 - L 2*2,- L 22|- L 2i2,- L 22 

p = -2 m * a -^92° --gsv -7, m v* a + 2 m ^o + 2 m ^° + 2 mpPo3 

i ^ 2 j b + 1 r kp B k 4 dk r Kp i k 4 dk 

+ 3 B ^ p ^^J {k 2 +m *Ky / 2 + 3 i= ^-*J (fc2 + m 2)l/2 • ( ) 

The mixed phase of hadronic and K~ condensed matter is governed by the Gibbs condi- 
tions for thermodynamic equilibrium along with global charge and baryon number conser 
vation laws 



35 



451 ]. The Gibbs phase rules read, 

P h = P K , (25) 

A = /4, (26) 

where fi% and /ijj are chemical potentials of baryons B in the pure hadronic and K~ con- 
densed phase, respectively. The global charge neutrality and baryon number conservation 
laws are 

(i - x)Q h + xQ K = 0, (27) 



n B = (1 - X )n h B + X">b , ( 28 ) 
where x is the volume fraction of K~ condensed phase in the mixed phase. The total energy 
density in the mixed phase is given by 

6 = (1 - X )e h + xe K ■ (29) 

Pressure and density variations associated with the r-mode oscillation drive the system 
out of chemical equilibrium. Microscopic reaction processes restore the equilibrium. Weak 
interaction processes are most important in this situation. The corresponding bulk viscosity 



coefficient cou 
lfj 
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d play a significant role on the mode. Earlier it was shown by various authors 



15l | that non-leptonic processes involving hyperons might lead to a high value 



for the bulk viscosity coefficient. The real part of 
terms of relaxation times of microscopic processes 



12 



k viscosity coefficient is calculated in 



46| 



1 + (W) 

where P is the pressure, r is the net microscopic relaxation time and 7^ and 70 are 'infinite' 
and 'zero' frequency adiabatic indices respectively. The factor 

P dn n dnb ' 

can be determined from the equation of state (EoS). Here x n = ^ gives the neutron fraction 
in the equilibrium state and = YIb 71 ^ is the total baryon density. In the co-rotating frame, 
the angular velocity (u) of (l,m) r-mode is related to angular velocity (Q) of a rotating 
neutron star as u = & [41] • 

The partial derivatives in both phases are calculated using the Gibbs-Duhem relation. In 
the pure hadronic phase, this relation gives P h = n 1 ^^ + vt^l + n e \i e + n^fi e — e h . So the 
partial derivative is given by 

dPk -u h + n^ + n^-^ (32) 



dn h n rn n dn h n p dn h n dn\ 



Using \i h n = f^, we obtain 



Qph 



()nh ri h y nn + n h p a h pn} (33) 

where aia is defined by a™ = ( ) . Similarly, in the pure condensed phase, the partial 
derivative is 

QpK 

-q^K = H n a nn + n p ®pn + ^K~0^- n ■ (34) 



In the mixed phase, we obtain the relation 

dP dP h 



dn r . 



dnl 



+ 



dP 



K 



(35) 



where n n = (1 - x)n h a + x n n ■ 

We are interested to study the bulk viscosity coefficient and corresponding damping 
timescale due to the non-leptonic weak decay process (21). For the reaction rate of this 
process, we may express all perturbed quantities in terms of the variation in neutron number 
density (n^) in the Bose-Einstein condensed phase. The relaxation time (r) for the process 
is given by 12] 

II' An 

(36) 



1 

T 



T K 8jj, 



Here, 5n^ = — is the departure of neutron fraction from its thermodynamic equi- 
librium value in the K~ condensed phase. The reaction rate per unit volume can be 
written as 



K 



(27T) 



d 3 ki d 3 k 2 d 3 k 3 



{\M\ 2 )5\h -k 2 - k 3 )F(E i )S°(E 1 -E 2 -E 3 



(37) 



](2nY J E 1 E 2 E 3 

where indices i = 1 — 3 refer to n, p and K~ respectively. Here we are considering K~ 
mesons in the s-wave (k 3 =0) condensate. So the conservation of momentum gives k\ = k 2 . 
Assuming that the matrix element has no angular dependence and performing the integrating 
over d 3 k 3 , we get 



K 



8(2vr) & E 3 



\k 1 \E 1 dE 1 \k 2 \E 2 dE 2 



E, 



{\M\ 2 )F{E i )8°{E 1 -E 2 - E 3 )dfl 1 dfl 2 



where F(Ei) is the Pauli blocking factor given by 

F(E i ) = f 1 (l-f 2 )-(l-f 1 )f 2 



(38) 



(39) 



and fi 



i 



. |A^| is the squared matrix element and (.) denotes the averaging over 

(l+e fcT 1 ) 

initial spins and sum over final spins. In this calculation fermions which reside close to their 
Fermi surfaces, contribute to the reaction. Integrating over E\ and substituting E2 ^ 2 = y 
and dE 2 = kTdy, we obtain 

• + °° dy 



r 



K 



(\M\ 2 ) 
16tt 3 £ 3 N 



ki 



Finally we arrive at the result, 



r 



K 



kT 



~ {\M\ 2 ) 2 
1Qit 3 E 3 n 



1 + e-y)(l + ev) 



5fi . 



(40) 



(41) 



8 



Here the in- medium energy of K~ mesons in the condensate is E% = 113 = fix- and k Fn is 
the Fermi momentum for neutrons in the condensed phase. It is noted that the reaction 
rate per unit volume is proportional to the overall chemical potential imbalance 5 fx = 8fx% — 
5fip — Sfj,K-- Now the relaxation time (36) can be written as 



1 

T 



(\M\ 2 ) 



5fx 
6^ 



16tt 3 ix k - ■ ' ' 12 ' 

Next we focus on the evaluation of the matrix element of the non-leptonic process (21). 
In the general form, this can be written as 

M=u(k 2 )(A + B l5 )u(k 1 ) , (43) 

where u{k{) and u{k 2 ) are the spinors of neutrons and protons respectively. Here A is the 
parity violating amplitude and B is the parity conserving amplitude. The squared matrix 
element is given by 

(44) 



(\M\ 2 ) = 2[(h ■ k 2 + mfmf) \A\ 2 + {k x ■ k 2 - mfmf) \B\ 2 ] . 

As fermion momenta lie close to their Fermi surfaces, we can write k\ ■ k 2 
\k± | \k 2 1 cos 9 = Li^lXp — kp n . This leads to the squared matrix element 

(\M 2 \) = 2[(ii*rf - k Fn k Fp + mfmf) \A\ 2 + (rfrf - k Fn k Fp - mfmf) \B\ 2 



E X E 2 — 



(45) 



Here we apply the weak SU(3) symmetry to the non-leptonic weak decay amplitudes for 
the process (21). Earlier weak decays of the octet hyperons were described by an effective 



SU(3) interaction with a parity violating (A) and parity conserving (B) amplitudes 49l. |50| 



The weak operator is proportional to Gell-Mann matrix A6 to ensure hypercharge violation 
| Ay | = 1 and |AJ| = 1/2 rule. The amplitudes were determined experimentally from 
the weak decay parameters of the octet hyperons 50|. We follow the same mechanism 
to determine the amplitudes in (21). The amplitudes are A = —1.62 x 10 -7 and B = 
-7.1 x 1(T 7 . 

The computation of the relaxation time is complete with the calculation of 4^. The 
chemical potential imbalance due to the non-leptonic process (21) is given by 



5 ix = 5[x n — 5ix p — 5ix K - 

= ( a nf U n + a np Sn p + U nK -8n K ~) 

- (®pJ n n + rfp^p + a pK-$n K -) 

- {oiK-nSn^ + a K - p 5nf + a K - K -5n K -) 



(46) 



We express S/j, in terms of 5n^ and obtain from the following constraints, 

n 

(1 - x)(K + 5n h p ) + X {K + 8rf) = , 
(l- x )5n h p + X (Sn« -6n K -) = , 

r h <; K 

6 fi p = dfi p , 

Sfi = Stf . (47) 

In the above constraints, we have 5x = because number densities deviate from their 
equilibrium values not by changing the volume of the fluid element but by internal reactions 



12j . First two constraints result from the conservation of baryon number and electric charge 
neutrality. The other constraints are due to the equality of neutron and proton chemical 
potentials in the hadronic and condensed phases. Further the last two constraints in (47) 
can be written as 

{a h p Jn h n + a h pp 6n h p ) - (a«Jn« + agSnf + a pK -5n K -) = , 
(a h nn 8n h n + a h np 5nl) - (a^Jn* + a« p 5n« + a nK -5n K -) = . (48) 

Next we calculate in the hadronic as well as K~ condensed phases using the EoS. 
These quantities in the hadronic phase are given by 

d» h n 



nn 



h (g^V 1 (9 p n\ 2 mfg aN da 



where 



da 



9aN 

mi 



dn h n D 



da 



9 a N _^ ====== 



dn h p D 



(51) 



(52) 



and 



g=l + ij^ + -£ gfe+il (05) 2 . (53) 
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Also, «p p and a pn are obtained by interchanging n ^ pin a\ n and a^ p respectively. Similarly, 
atij in the K~ condensed phase are given by 

k _ fyfi ( 9^n\ 2 1 / g P N\ 2 7T 2 m*?g aN da 



k_(9un\ 2 1 (9pn\ 2 rn*tfg aN da 
<9)"n / 9^n9^k\ . 1 / g P Ng P K\ m*fg aN da 

°- = = - 1— J + 3 1— J - + mjr . • (56) 

fyfg^ __ / 1 / g P Ng P K \ da 



^ _ / nu)i\nu)i\ J _|_ _ i ypjvi/pj\ i ^ (57) 

<9n^ \ m 2 / A\m 2 p J ° dn% 

2 / \ 2 / \2 

ffyiR-- / \ , / g<pK \ . 1 / g P K \ da da* 



OLK-K- = ~ = + +T - 9aK-R 9a*Kn ' ( 58 ) 

dn K - \m u J V m / 4 \ m p/ dn K - dn K - 



where, 



, * K 



do- ( g aN \ V fc 



dn^ \ mi ) D 

o 9<jK 

oa m 2 



■ 



dn K - D' 
da* g a * K 
dn K - m 2 * ' 



(59) 



and J^jr = = 0. Here D' has the same form as D in (53), but D' has to be calculated 
using the meson fields in the condensed phase. The coefficients , cxjf n , cupx- and a K - p 
are obtained by interchanging n <-> p in , a^ p , and respectively. 

The bulk viscosity damping timescale (r B ) due to the non-leptonic process (21) is defined 

by 

1 1 dE 

(60) 

where E is the energy of the perturbation in the co-rotating frame of the fluid and is given 



by 



12| 



1 r R 

E = -a 2 Q 2 R~ 2 / e{r)r 6 dr. (61) 

Here a is the dimensionless amplitude of the r-mode, R is the radius of the star and e(r) is 
the energy density profile. The rate of energy dissipation is given by 

*E = -4tt [ R C(r) < | V • Sv\ 2 > r 2 dr, (62) 



dt 



o 



11 



where the angle average of the square of the hydro dynamic expansion 15|, |5jj is (|V • 
5v\ 2 ) = ^*g9Q (^) 6 1 + 0.86 (^) 2 (jtu^J an d ^ i s the mean energy density of a non-rotating 
star. Equation (60) is the bulk viscosity contribution to the imaginary part of the r-mode 
frequency. We also take into account time scales associated with gravitational radiation 
{t~gr), bulk viscosity due to modified Urea process (jy) involving only nucleons and the 
shear viscosity (tsv) and define the overall r-mode time scale (r r ) as 



1 1111 
— = + — + — + 

T~r T GR T B Tjj T SV 



(63) 



The gravitational radiation timescale is given by 



17, 



131072tt 



Q 6 



R 



e(r)r 6 dr . 



(64) 



t gr 164025 

The time scale corresponding to the bulk viscosity due to the modified Urea process (tjj) 
involving only nucleons is calculated in the same footing as the bulk viscosity damping time 



scale due to the non-leptonic process but using the fol 
coefficient for the modified Urea process given by 3], 



owing expression for the bulk viscosity 
53 



(u = 6 x 10 



25 



10 15 g/cm 3 



[WkJ (ilk) Zg,cms 



(65) 



The damping time scale due to the shear viscosity is given by 



,7] 



T~SV 



rj(r)r 4 dr 



R 



e(r)r G dr 



(66) 



Different scattering processes contribute to the shear viscosity (77) depending on whether 
neutrons and protons are superfluid or not ^, [l8]. In this calculation, neutrons and pro- 
tons are both normal fluids. The most important contribution comes from neutron-neutron 
scattering in this case and the shear viscosity coefficient has the form 



V 



2 x 10 18 



9/4 



\10 9 K J 



g/cms . 



(67) 



10 15 5'/cm 3 / 

The total r-mode time scale is a function of angular velocity, neutron star mass and 
temperature. Therefore, solving y = 0, we calculate the critical angular velocity above 
which the r-mode is unstable whereas it is stable below the critical angular velocity. 
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III. RESULTS AND DISCUSSION 



We need to know meson-nucleon and meson-kaon coupling constants for this calculation. 
The nucleon-meson coupling constants are determined from the nuclear matter saturation 
properties of binding energy E/B = —16.3 MeV, baryon density n = 0.153 fm -3 , asym- 
metry energy coefficient a asy = 32.5 MeV, incompressibility K = 240 MeV and effective 



nucleon mass m* N /m^ = 0.70. The coupling constants are taken from Ref. 53J and this set 
is known as GM. We determine kaon-meson coupling constants using the quark model and 
isospin counting rule. The vector coupling constants are given by 

9uK = t^QuN and g pK = g pN . (68) 

The scalar coupling constant is obtained from the real part of K~ optical potential depth 
at normal nuclear matter density 

Ur {no) = -g<rK(y - g^K^o ■ (69) 

It was found in various calculations that antikaons experienced an attractive potential and 



54 



55 



56 



57 



58, 



591 ] . The strength 



kaons had a repulsive interaction in nuclear matter 
of antikaon optical potential ranges from shallow (—40 MeV) to strongly attractive (—180 
MeV). Moreover, it was noted that K~ condensation could be a possibility for antikaon 
potential ~ 100 MeV or more attractive. Here we perform the calculation for a set of values 
of antikaon optical potential depths at normal nuclear matter density. We obtain kaon- 
scalar meson coupling constants g a x = 2.1914,3.0426,3.8937 corresponding to Ur{tiq) = 
-100,-120,-140 MeV. 

The strange meson fields couple with (anti)kaons. The a*-K coupling constant is 
g a *K — 2.65 as determined from the decay of /o(925) and the vector <p meson coupling 



with (anti)kaons V^g^K — 6.04 follows from the SU(3) relation 

We begin our investigation with K~ condensation in nuclear matter and denote this sys- 
tem as iapK~ . Baryon densities corresponding to lower (ui = rii/n ) and upper (u u = n u /n ) 
boundaries for K~ condensation with antikaon optical potential depths at normal nuclear 
matter density Ux(n ) = —100, —120, —140 MeV are shown in Table I. The K~ condensa- 
tion sets in at a density ni = uiUq and a pure K~ condensed phase with nucleons embedded 
in it begins at a density n u = u u n . It is evident from the Table that there is no mixed 
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phase for U^(jio) = — 100 MeV. We further note that the phase transition is of second order 
for \Uj^{no)\ < 100 MeV. On the other hand, the onset of K~ condensation shifts towards 
lower density as the antikaon optical potential becomes more attractive. The mixed phase 
is found to be wider for Uj^(n ) = —140 MeV. These findings were already reported in Ref. 
35( | . Another interesting problem arises when hyperons a ppe ar around threshold densities 



of K~ condensation. We showed in an earlier calculation [16| by extending the Lagrangian 
density in (1) to include hyperons that A hyperons first appeared around baryon density 
2.58n . However, X hyperons are excluded from the system because of a strongly repulsive 
E-nuclear matter interaction whereas massive H hyperons appear at higher densities. In 
this context, we continue our calculation with K~ condensation in npA matter denoted by 
npAA' - . In this case, A-meson couplings are taken from Ref. [3]. The threshold density 
(Ut h ) for the appearance of A hyperons are listed in Table I. For C/#-(no) = —100, —120 MeV, 
A hyperons appear before the K~ condensate. Consequently, K~ condensation is delayed 
to higher densities. Those values are recorded in parentheses. On the other hand, the early 
occurrence of K~ condensation for [/g-(no) = —140 MeV delays the appearance of A hyper- 
ons to higher density. Therefore, the composition of matter and the strength of antikaon 
optical potential depth determine the onset of K~ condensate and A hyperons in npAK~ 
matter. It was already noted that the bulk viscosity coefficient due to the non-leptonic 
process involving A hyperons had a very large value. In this calculation, we study the bulk 
viscosity due to non-leptonic process (21) in npfC~ matter and find out whether kaon bulk 
viscosity alone can damp the r-mode without hyperons. 

Now we describe the composition of /5-equilibrated and charge neutral matter in the 
presence of K~ condensation for £7^-(n ) = — 120 MeV. Populations of various particles are 
displayed in Fig. 1. Before the appearance of K~ condensate, particle fractions increase with 
baryon density. In this situation, charge neutrality is maintained among protons, electrons 
and muons. The mixed phase begins with the onset of K~ condensation at 3.26 n Q . As 
soon as K~ condensate is formed, its rapid growth replaces electrons and muons from the 
system. Consequently, the proton density becomes equal to the density of K~ mesons in 
the condensate. 

The equations of state, pressure versus energy density, with and without K~ conden- 
sate are shown in Fig. 2. The solid line represents the overall EoS without K~ condensate 
whereas the long dashed, short dashed and dotted lines correspond to those with the conden- 
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sate for Ug = —100, —120, —140 MeV respectively. We have already noted that the phase 
transition is of second order for U^{no) = —100 MeV. On the other hand, we find that two 
kinks on the EoS involving the condensate mark the beginning and end of the mixed phase 
for Uk(tio) = —120, —140 MeV. The EoS becomes softer in the presence of K~ condensate 
compared with the case without the condensate. The deeper the antikaon optical potential 
depth, the softer is the EoS. 



Next we construct models of non-rotating and uniformly rotating neutron stars [60j using 
EoS as shown in Fig. 2. The gravitational mass for non-rotating neutron star sequence 
and the sequence of neutron stars rotating at their Kepler frequencies are plotted with 
central energy density e c in Fig. 3. The inclusion of K~ condensate in EoS for different 
values of U^(uq) reduces maximum masses for non-rotating as well as rotating neutron 
stars. The maximum gravitational masses, equatorial radii, corresponding central baryon 
and energy densities and Kepler periods for the mass shedding limit are given by Table 
II. The corresponding values for non-rotating neutron stars are shown in parentheses. It 
is evident from Table II that the softest EoS corresponding to U^{no) = —140 MeV leads 
to a maximum non-rotating neutron star mass 1.343 M which is well below the most 
accurately measured Hulse- Taylor pulsar mass 1.44 M & . Therefore, this EoS is ruled out by 
observations. We adopt antikaon optical potential depths U^no) = —100, —120 MeV for 
the calculation of bulk viscosity. 

The equation of state enters as an input in the calculations of adiabatic indices, relaxation 
time and bulk viscosity coefficient. In particular, partial derivatives of pressure and chemical 
potentials with respect to neutron number density are calculated using the EoS. Figure 4 
exhibits the difference of fast and slow adiabatic indices as a function of normalised baryon 
density for U^(no) = —100, —120 MeV. We find jumps in the difference of adiabatic indices 
(7oo — 7o)- For Ux(n Q ) = —100 MeV, this jump occurs at the onset of K~ condensation 
whereas these are attributed to kinks on the EoS for Ug = —120 MeV in Fig. 2. Those kinks 
lead to discontinuities in the partial derivatives in Eq. (31). This kind of jump in adiabatic 
indices was earlier noted in quark-hadron phase transition in neutron star matter by others 



6l|. We further note that the value of (700 — 70) is suppressed for Ux(n ) = —100 MeV 



compared with the case of Uj^(n Q ) = —120 MeV. The relaxation time for the non-leptonic 
weak process (21) is plotted with normalised baryon density in Fig. 5 for different values 
of antikaon optical potential depths. In this case, the maximum value of relaxation time 
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is ~ 10 _11 s for Ux{no) = —120 MeV. This value is a few order of magnitude smaller than 
that of the non-leptonic process involving non-superfluid A hyperons jig, Also, we note 
that the relaxation time given by Eq. (42) is independent of temperature because there is 
no temperature dependence in the reaction rate in (41). This may be attributed to the fact 
that the K~ condensate is treated here as a zero temperature Bose-Einstein condensate. On 



the other hand, the relaxation time for the non 
was found to be inversely proportional to T 2 12 



eptonic weak process involving hyperons 
I5I . It is worth mentioning here that 



the role of thermal kaon processes in the non-equilibrium dynamics of kaon condensation 
was studied by others [62j, |63J] . 

The bulk viscosity coefficient in the kaon condensed phase is shown as a function of 
normalised baryon density in Fig. 6. The factor ujt in the bulk viscosity coefficient given 
by Eq. (30) is negligible compared with unity over the whole range of baryon densities 
considered here. We note that the bulk viscosity coefficient is discontinuous at the phase 
boundaries of hadronic and K~ condensed matter. This is attributed to kinks in the EoS and 
discontinuities in (7^ — 70) . It is worth mentioning here that the bulk viscosity coefficient 
in the kaon condensed matter is several orders of magnitude smaller than the non-superfluid 
hyperon bulk viscosity coefficient [l2,[l^, We also note that the bulk viscosity coefficient 
is larger for U^(no) = —120 MeV than that of the case with £7^-(no) = —100 MeV in the 
lower density regime. We use the bulk viscosity coefficient corresponding to Ux(no) = —120 
MeV for the calculation of damping time scale and critical angular velocity. 

The computation of damping time scale due to bulk viscosity given by Eq. (60) involves 
the energy density profile, kaon bulk viscosity profile and structures of rotating neutron 
stars. The sequence of neutron stars rotating at Kepler frequencies is calculated by using 
the model of Stergioulas 60]. We have already shown the mass shedding limit sequences of 
rotating neutron stars in Fig. 3. For the calculation of damping time scale, rotating neutron 
stars which contain K~ condensate in its interior are to be considered. The appearance of 
K~ condensate is sensitive to the central density of a rotating neutron star which ,in turn, 
depends on the rotation period of the star. It is worth mentioning here that the canonical 
neutron star of 1.4M rotating at the Kepler period 1.22 ms has a central baryon density 
2.20n and that of its non-rotating counterpart has a central baryon density 2.76n for 
U^in^) = —100, —120 MeV. This shows that central densities of 1.4M canonical neutron 
star are well below threshold densities of K~ condensation. To demonstrate the effect of 
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kaon bulk viscosity on the r-mode instability, we choose a neutron star of gravitational 
mass 1.63M having central baryon density 3.94 no and rotating at an angular velocity 
Vl rot = 1180s -1 from the sequence of rotating neutron stars corresponding to Ug = —120 
MeV. This neutron star contains K~ condensate in its core because the central baryon 
density is well above the threshold density of K~ condensation. The energy density profile 
of this rotating star is plotted as a function of equatorial distance in Fig. 7. Similarly, the 
kaon bulk viscosity profile of this neutron star as a function of equatorial distance is shown 
in Fig. 8. Here we note that the bulk viscosity profile drops to zero value beyond 2 km. This 
happens because the baryon density beyond 2 km decreases below the threshold density of 
K~ condensation and the non-leptonic process (21) ceases to occur there. 

Using the energy density and bulk viscosity profiles, we estimate the damping time scale 
corresponding to the bulk viscosity due to the non-leptonic process (21). We also consider 
the contributions of modified Urea process involving nucleons as well as the shear viscosity 
to the total r-mode time scale as given by Eq. (63). We obtain critical angular velocities 
as a function of temperature solving l/r r = for a rotating neutron star mass 1.63M Q 
corresponding to the case Uk{tiq) = —120 MeV and show this in Fig. 9. The r-mode is 
unstable above the critical angular velocity curve and stable below it. The bulk viscosity due 
to the modified Urea process dominates at higher temperatures and might damp the r-mode 
instability. The long dashed line which merges with the solid line above 5xl0 9 K, shows the 
modified Urea contribution balancing the gravitational radiation reaction. As temperature 
decreases, the bulk viscosity due to the modified Urea process also decreases. Consequently 
the corresponding damping time scale due to the modified Urea process is longer than the 
gravitational radiation growth time scale and the r-mode instability is not suppressed below 
10 9 K by this process. Now we consider the shear viscosity. As the shear viscosity coefficient 
is proportional to T -2 , it might suppress the instability at lower temperatures. This is shown 
by the short dashed line which merges with the solid line below temperature 5xl0 9 K. Next 
we discuss what is the influence of kaon bulk viscosity on the r-mode instability. We know 
already that kaon bulk viscosity is independent of temperature. The bulk viscosity damping 
time scale (rg) due to the non-leptonic process (21) and the time scale associated with 
gravitational radiation are also independent of temperature but depend on angular velocity 
(Q). Further investigation reveals that the damping time scale tb is always longer than the 
gravitational radiation growth time scale tqr for all values of Q ranging from to Vt rot . This 



17 



implies that the bulk viscosity due to the non-leptonic process in K~ condensed matter can 
not damp the r-mode instability. When all viscous processes are taken into account, we find 
an intermediate temperature regime where neither shear viscosity nor bulk viscosity can 
damp the r-mode instability, is marked by the solid line. The critical angular velocity curve 
is determined due to the interplay of tjj, t sv and t gr . 



IV. SUMMARY AND CONCLUSIONS 



We have investigated the role of K~ condensation on bulk viscosity and r-mode instability 
in neutron stars. The bulk viscosity coefficient and the corresponding damping time scale 
due to the non-leptonic process n # p + K~ have been calculated. In this calculation, we 
have considered first and second order kaon condensation and constructed equations of state 
for different values of antikaon optical potential depths within the framework of relativistic 
field theoretical models. It is noted that kaon bulk viscosity and the corresponding damping 
time scale are independent of temperature. We find that the bulk viscosity coefficient in 
K~ condensed matter is suppressed by several orders of magnitude compared with the non- 
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161 ] . Further we note that kaon bulk 



superfluid hyperon bulk viscosity coefficient] 
viscosity due to the process (21) alone can not damp the r-mode instability because kaon 
bulk viscosity damping time scale is always longer than the gravitational radiation growth 
time scale. 

We have also discussed the role of hyperons on the appearance of K~ condensate. The 
early appearance of A hyperons delays the onset of K~ condensation to higher densities for 
antikaon optical potential depths U R (n ) = -100, -120 MeV. For U R (n ) = -140 MeV, K~ 
condensation occurs before A hyperons are populated in the system. The latter situation 
was already reported in Ref. 42| for a different parameter set. In this connection, we note 
that the large value of hyperon bulk viscosity might damp the r-modes in neutron stars 
when hyperons appear in neutron stars. On the other hand, the r-modes in neutron stars 
containing only kaon condensate could be damped by the bulk viscosity of modified Urea 
process at higher temperatures and by the shear viscosity at lower temperatures as it is 
demonstrated in Fig. 9. It is to be noted here that the onset of K~ condensation is very 
sensitive to the rotation period of the neutron star under investigation. The gravitational 
instability window for the case with K~ condensate in nuclear matter would be different 
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from the situation when hyperons are included with the condensate. 
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TABLE I: For neutron star matter with nucleons and antikaon condensate (npK~), baryon den- 
sities u = n/no, corresponding to the lower (v!) and upper (u u ) boundaries of the mixed phase 
for first order K~ condensation at various values of antikaon optical potential depths £7^ (no) (in 
MeV) at the saturation density no = 0.153 fm~ 3 are given by this table. Threshold densities of 
A hyperons, u£ h , are shown for neutron star matter composed of nucleons, A hyperons and K~ 
condensate. The values in parentheses are obtained in the presence of A hyperons. 





u l 


u u 


< 


-100 


3.88 (4.70) 


3.88 (4.70) 


2.58 


-120 


3.26 (3.47) 


4.62 (5.35) 


2.58 


-140 


2.46 (2.46) 


5.40 (6.07) 


3.91 
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TABLE II: Maximum masses M max , equatorial radii R and their corresponding central baryon 
densities ti CC nt = n cent/^o an d energy densities e c for nucleons-only (np) star and for stars with 
further inclusion of antikaons (npK~) rotating at Kepler periods Pr are listed below. Results at 
different values of antikaon optical potential depths [7^ (no) (in MeV) at the saturation density are 
recorded here. The values in parentheses correspond to non-rotating neutron stars. 



EoS 




P K (ms) 


^cent 


£ c (10 15 g/cm 3 ) 


M max /M 


R(km) 


np 




0.66 


6.40 (7.00) 


2.15 (2.43) 


2.371 (2.019) 


14.84 (11.02) 


npK~ 


-100 


0.80 


5.33 (6.12) 


1.66 (1.96) 


2.189 (1.817) 


16.67 (11.98) 




-120 


0.91 


4.62 (5.20) 


1.38 (1.58) 


2.005 (1.642) 


17.73 (12.62) 




-140 


1.06 


3.99 (4.56) 


1.15 (1.33) 


1.648 (1.343) 


18.51 (12.98) 



23 




Fig. 1. Particles abundances are plotted with normalised baryon density for antikaon optical 
potential depths at normal nuclear matter density U^tiq) = — 120 MeV. 
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Fig. 2. The equation of state, pressure P vs energy density e , for nucleons-only matter with and 
without K~ condensate. The results including the condensate are for antikaon optical potential 
depths at normal nuclear matter density U^(uq) = —100, —120, —140 MeV. 
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Fig. 3. The gravitational mass for non-rotating sequence as well as mass shedding limit sequence 
of rotating neutron stars are plotted with central energy density for the equations of state shown 
in Fig. 2. Different lines have same meaning as in Fig. 2. 
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Fig. 4. The difference of adiabatic indices, 7^ — 70, is shown as a function of normalised baryon 
density for EoS including K~ condensate with antikaon optical potential depths at normal nuclear 
matter density U R (n ) = -100, -120 MeV. 
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Fig. 5. Relaxation time is plotted with normalised baryon density for the non-leptonic process 
(21) and antikaon optical potential depths at normal nuclear matter density U^(no) = —100, —120 
MeV. 
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Fig. 6. Bulk viscosity coefficient is exhibited as a function of normalised baryon density for the 
process (21) and antikaon optical potential depths at normal nuclear matter density U^uq) = 
-100,-120 MeV. 
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Fig. 7. Energy density profile is shown with equatorial distance for a rotating neutron star of 
mass 1.63 M Q corresponding to the EoS with antikaon optical potential depth at normal nuclear 
matter density U^(n ) = -120 MeV. 
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Fig. 8. Bulk viscosity profile is shown with equatorial distance for a rotating neutron star of mass 
1.63 M Q . The antikaon optical potential depth is same as in Fig. 7. 
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Fig. 9. Critical angular velocity for 1.63 M neutron star is plotted as a function of temperature. 
The antikaon optical potential depth is same as in Fig. 7. The solid line denotes the angular 
velocity curve when all viscous processes are included. The contributions from the modified Urea 
bulk viscosity and shear viscosity are shown by long and short dashed lines which merge with the 
solid line at higher and lower temperatures respectively. 
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